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ABSTRACT 

Context. GJ1214b, the 6.55 Earth-mass transiting planet recently discovered by the MEarth team, has a mean density of ~35% of that 
of the Earth. It is thought that this planet is either a mini-Neptune, consisting of a rocky core with a thick, hydrogen-rich atmosphere, 
or a planet with a composition dominated by water 

Aims. In the case of a hydrogen-rich atmosphere, molecular absorption and scattering processes may result in detectable radius 
variations as a function of wavelength. The aim of this paper is to measure these variations. 

Methods. We have obtained observations of the transit of GJ12I4b in the r- and I-band with the Isaac Newton Telescope (INT), in 
the g-, r-, i- and z-bands with the 2.2 meter MPI/ESO telescope, in the Kj-band with the Nordic Optical Telescope (NOT), and in the 
Kr-band with the William Herschel Telescope (WHT). By comparing the transit depth between the the different bands, which is a 
measure for the planet-to-star size ratio, the atmosphere is investigated. 

Results. We do not detect clearly significant variations in the planet-to-star size ratio as function of wavelength. Although the ratio at 
the shortest measured wavelength, in g-band, is 2tr larger than in the other bands. The uncertainties in the K, and bands are large, 
due to systematic features in the light curves. 

Conclusions. The tentative increase in the planet-to-star size ratio at the shortest wavelength could be a sign of an increase in the 
effective planet-size due to Rayleigh scattering, which would require GJ1214b to have a hydrogen-rich atmosphere. If true, then the 
atmosphere has to have both clouds, to suppress planet-size variations at red optical wavelengths, as well as a sub-solar metallicity, 
to suppress strong molecular features in the near- and mid-infrared. However, star spots, which are known to be present on the host- 
star's surface, can (partly) cancel out the expected variations in planet-to-star size ratio, because the lower surface temperature of the 
spots causes the effective size of the star to vary with wavelength. A hypothetical spot-fraction of ~10%, corresponding to an average 
stellar dimming of ~5% in the i-band, would be able to raise the near- and mid-infrared points sufficiently with respect to the optical 
measurements to be inconsistent with a water-dominated atmosphere. Modulation of the spot fraction due to the stellar rotation would 
in such case cause the observed flux variations of GJ 12 14. 

Key words, techniques: photometric - stars: individual: GJ1214 - planetary systems 



1. Introduction 

During the transit of an exoplanet, the light from its host-star 
filters through the planet's atmosphere, and atmospheric signa- 
tures from molecules and atoms get imprinted on the transmis- 
sion spectrum. The strength of the atmospheric features is de- 
pendent on the atmospheric scale height H, which in turn de- 
pends on the temperature T, the mean molecular weight fi and 
the planet's surface gravity g as H=^. For hot-Jupiters this 
scale height is a few hundred kilometers. For instance, the ex- 
pected atmospheric scale height of HD 189733b for an H2 dom- 
inated atmosphere is ~200 km, while for HD209458b this is 
about 700 km. For transit observations, the observable param- 
eter we are interested in is ARp/R, ~H/R,, which is ~0.0004 for 
HD 189733b, and -0.0009 for HD209458b. Nevertheless, signa- 
tures from the atmospheres of hot-Jupiters have been detected 



* Photometric timeseries are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
[http://cdsw eb.u-strasbg.fr/cgi-bin/qc at ?J/A+A/| 



in transmission spectra, bo t h from atoms, includ i ng sodium (eg. 
'Charbonneau et al." '2002'; Snellen et al. 2008; Redfiel d et al] 
2008) . potassium (Sing et al. 2()lla; Colo net al. 201( |, hy- 
droge n (IVidal-Madiar et al.l I2OO3I). "carbon ([V idal-Madi ar et all 
2004') and oxvgen (IVidal-M adiar et al."2004'). as well as from 
molecules such as water (T inetti et al. 2007), althou gh this 
is disputed by iDesert et~al.l (120091) . methane (ISwain et al.l 
2008.). which has been challenged by ISing et al.1 (l2009h an d 
iOibson etal.l(l20TTl) . and carbon monoxide ([Snellen et alj2010l) . 
In addition, a gradual increase of the planet-to-star radius ra- 
tio of HD 189733b has been detected toward shorter wave- 
length vwhkhhas_beenjt^^ by haze par- 
ticles (iPont et al.ll20"08llSing et alJl201 Ibh . 

For cooler and smaller planets, the scale height decreases, 
and the AR^/R* becomes much smaller. For the Earth the scale 
height is only 8.5 km, which corresponds to a change in the 
radius ratio of ~10"^, which is very small. However, f or the 
recently discovered super-Earth GJ1214b (Charbonnea u et al.[ 
I2009i) . the scale height can be similar to that of HD 189733b, 
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due to the low (~Q.35pearth) density of the planet, if the atmo- 
sphere is dominated by H2. Since its host-star, GJ1214, is ~4 
times smaller than the host-star of HD 189733b, the change in 
the planet-to-star radius ratio is ~4x larger at ARp/R, =0.0016. 
This makes GJ1214b an ideal candidate to search for the signa- 
tur es of its atmosphere, desp ite its small mass and radius. 

[Rogers & Seager (2010) presented three formation scenar- 
ios for GJ1214b that explain its low density. These scenar- 
ios also provide predictions on the composition of the plan- 
ets atmosphere. If GJ1214b formed as a predominantly water- 
rich planet, the atmospheric scale height would be small, since 
the mean molecular weight of water is high. However, if 
GJ1214b's low density is due to out-gassing from a rocky 
planet, or due to it being formed as a mini-Neptune, the atmo- 
sphere is expected to consist predominantly of hydrogen and he- 
liu m, with a low mean mo l ecular weight. Atmospheric models 
by Miller-Ri cci & FortnevI ( 1201 Ol) showed that it is possible to 
get detectable signatures from an atmosphere with a large scale 
height, especially from molecules in near-infrared, but also from 
the scattering of light in the optical part of the spectrum. 

Recently, the first transit transmission results for the at- 
mos phere of G J1214b have been presented in the litera- 
ture. iBean et al] (12010 ) found that their spectroscopy in the z- 
band showed no sign for the presence of a thick, hydrogen-rich 
atmosphere, which argues for a water-rich atmosphere, some- 
th ing which i s also consistent with the mid-infrared observations 
of iDesert et al . (2011a) . These obs ervations are in contrast with 
the results from C roU et al.l who show that the transit in 

the Kj-band is deeper than the transit in the J-band, and therefore 
consistent with an atmosphere with a large scale-height and low 
mean-molecular weight. 

In this paper we present the results for our multiband transit 
photometry of GJ1214b, covering a wavelength range from the 
g-band in the optical to the Kc-band in the near-infrared. In sec- 
tion |2] we present our observations, followed by the data reduc- 
tion in section[3]and transit fitting in section|4] Subsequently we 
discuss the influence of stellar variability in section|5]and present 
and discuss the transmission spectrum of GJ1214b in section|6] 
Finally we give the conclusions in section]?] 
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Fig. 1. Raw light curves for the INT r-band (top panel) and the 
NOT Kj-band observations (bottom panel). In the bottom panel 
we have also overplotted the flux in the aperture due to the sky- 
background 



order to keep the peak count-levels in the linear regime of the 
detector. This has the added benefit that the light is spread over 
more pixels, reducing the impact of flat-fielding errors. 



2. Observations 

2.1. WFC observations 

A transit of GJ1214b was observed in Sloan r-band {Ac-621 nm) 
on May 26, 2010 between 2:57 UT and 5:15 UT with the Wide 
Field Camera (WFC) on the 2.5 meter Isaac Newton Telescope 
(INT). An exposure time of 60 seconds was used resulting in 
89 frames with an average cycle time of 93 seconds. Only the 
central detector of the WFC (CCD4) was used for the analysis, 
with a pixel scale of 0.33"/pixel, this CCD has a field of view of 
675" by 1350", sufficient to observe a number of reference stars 
of similar brightness simultaneously with the target. 

The moon was almost full and the conditions were strongly 
non-photometric, with the transparency dropping to below 20% 
during several frames (see top panel of Fig. [T]i. 

On July 29, 2010, a transit of GJ1214 was observed in I-band 
(/lc=822 nm) with the same instrument. The observations started 
at 21:23 UT and lasted for just over 3 hours. An exposure time 
of 50 seconds was used, resulting in a total of 142 frames with 
an average cycle time of 81 seconds. In this case the night was 
photometric. 

Since GJ1214b is about 5.5 times brighter in I-band than in 
r-band, we significantly defocused the telescope in the I-band in 



2.2. GROND griz-band observations 

On the night of July 3, 2010, we obtained simultaneous obser- 
vations of GJ1214b in the g {Ac-459 nm), r (Ac-622 nm), i 
iAc^l64 nm), and z-band (A r=^99 nm) with the GROND in- 
strument dGreiner et al.ll2008h on the 2.2 meter MPI/ESO tele- 
scope at La Silla in Chile. The field of view in each of the 
wavelength channels is 5.4' by 5.4', which is sufficient to ob- 
serve both GJ1214 and a set of reference stars simultaneously. 
The observations started at 00:16 UT and lasted until 04:06 UT. 
During this time we obtained 280 frames in each of the four op- 
tical bands. The exposure time was varied from 20 to 30 seconds 
to avoid saturation of the CCDs. The average cycle time was 50 
seconds. 

2.3. NOTCam Krband observations 

We obtained a band transit observations with the NOTCam 
instrument on the Nordic Optical Telescope (NOT) simultane- 
ously with our INT r-band observations on May 26, 2010. The 
observations were caiiied out in service mode, and the wide field 
imaging optics and a K^-band filter {Ac-2.15 fim) were used. 
The pixelscale of this setup is 0.234"/pixel, resulting in a field 
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of view of the detector of 4 by 4 arcminutes. This field of view 
is sufficient to allow simultaneous observations of one reference 
star of similar brightness to GJ1214 as well as a reference star 
that is 4x fainter than GJ1214. The field of view of the detector 
was rotated to make sure that bad regions on the detector were 
avoided for all three stars. Since the NOT is located on the same 
mountain as the INT, these observations suffer from the same 
non-photometric conditions, with the transparency dropping to 
25% for parts of the light curve (see Fig. [1). This strongly af- 
fects the observations, since the sky background dominates over 
the object flux for the larger apertures, especially during times 
of low transparency. 

Since GJ1214 is bright at near-infrared wavelengths, we de- 
focused the telescope in order to allow for the relatively long in- 
tegration time and reducing the sensitivity to flat-fielding errors, 
although this also increases the impact of the sky background. 

The exposure time was set to 4 seconds, to allow for rel- 
atively efficient observations, with the large overheads induced 
by the NOTCam system. The average cycle time was 16 seconds, 
allowing us to capture 518 frames in 2 hours and 25 minutes of 
observations. 

In order to increase the stability of the system for the obser- 
vations, and to decrease the telescope overheads, we observed in 
staring mode, with guiding, keeping the centroid of the star con- 
stant to within 4 pixels during the observations. Since this ob- 
servation strategy does not allow us to subtract the background 
from the images, we obtained a set of dithered observations after 
our transit observation, from which we constructed a background 
map. 

2.4. LIRIS Kc-band observations 

The Kc-band {Ac-2.21fim) observations were obtained 
with the Long-slit Intermedi ate Re s o lution Infrared 
Spectrograph (LIRIS; lAcosta-Pulido et al.1 l2002h on the 
William Herschel Telescope (WHT) on the night of August 
17, 2010. The pixel-scale of the LIRIS detector is 0.25"/pixel, 
resulting in a field of view of the detector of 4.2 by 4.2 arcmin- 
utes. The same reference stars as for the NOTCam observations 
were used. 

A similar observing strategy as that of the NOTCam obser- 
vations was followed, defocusing the telescope and performing 
the observations in staring mode. An exposure time of 10 sec- 
onds was used, resulting in an average cycle time of 12.5 sec- 
onds. We obtained 1010 frames during the 3 hours and 26 min- 
utes of observations, of which 300 were taken during the tran- 
sit. From these we discarded the first 40 frames, 28 of which 
were taken with the stars on a different position on the detector, 
and the others because they were taken with a different defo- 
cus of the telescope. The frames were obtained in sequences of 
100, from which the first 3 frames clearly suffered from the reset 
anomaly. These frames were also excluded from further analy- 
sis (27 frames in total). As for the NOTCam observations, we 
observed a field offset from GJ1214 in order to construct a back- 
ground map. 



3. Data reduction 



3.1. Optical data 

A standard data reduction was performed for the WFC-INT data: 
after bias subtraction the images were flatfielded with a twilight 
flat. In addition a fringemap was used to coiTect the WFC I-band 
data for strong fringing. 

Subsequently aperture photometry was performed on 
GJ1214 and a set of reference stars (24 for the r-band and 14 
for the I-band), using apertures of 10 and 19 pixels for the r- 
and I-band respectively. The background was measured in an- 
nuli centered on the stars, using an inner radius of 25 and an 
outer radius of 50 pixels for the WFC r-band and 35 and 70 pix- 
els for the I-band. 

Since the weather conditions during the INT r-band observa- 
tions were strongly non-photo-metric, extinction gradients could 
be present in the photometry of the stars across the detector. 
Rather than directly combining light curves of all stars to ob- 
tain a reference light curve, we fitted the normalised photome- 
try of the reference stars in each frame with linear function in 
the (x,y)-position and background flux. We then calculated the 
value of the best fit function for GJ1214, to remove the extinc- 
tion effects of the atmosphere. Subsequently we fitted an airmass 
curve to the corrected light curve of GJ1214 to adjust for second 
order colour effects. The resultant light curve is shown in the left 
panel of Fig. |2l The out-of-transit noise-level is 1.610 '' per data 
point, which is 1.3 times larger than expected from Poisson and 
read-noise statistics. 

Since the WFC I-band data were taken under photomet- 
ric conditions, we simply combined the photometry of all the 
reference stars for the reference light curve. Subsequently we 
corrected for the airmass as for the r-band data. The resultant 
light curve is shown in the right panel of Fig. |2] The out-of- 
transit noise-level is 4.5-10"^, while the expected noise level 
from Poisson and read-noise statistics is 3-10"^. However, this 
estimate does not take scintillation noise into account, which is 
difficult to quantify, and could easily dominate the noise at this 
level. 

The GROND data, taken in the g-, r-, i- and z-bands simulta- 
neously, were reduced in a similar way as the INT data. An aper- 
ture of 1 8 pixels was used for the aperture photometry, while the 
background was determined in an annulus around the star with 
an inner radius of 20 pixels and an outer radius of 30 pixels for 
all bands. 

Subsequently, a linear base-line was fitted to the light curves 
of each band resulting in an out-of-transit RMS of the light 
curves of 3.3■10"^ 1.6-10"^ 1.4-10"^ and 1.3- 10"^ for the g-,r-, 
i- and z-bands respectively. These noise-levels are significantly 
higher than expected from Poisson and read-noise statistics, pos- 
sibly due to the fact that these in-focus observations are much 
more sensitive to flatfielding inaccuracies. 

3.2. Near-infrared data 

As a first step in the NOTCam Kj-band data reduction a non- 
linearity correction was made to all frames on a pixel-by-pixel 
basis, using a non-linearity measurement from the NOTCam cal- 
ibration archive Q. The science frames were subsequently flat- 
fielded using a twilight flat. Those frames obtained for sky- 
subtraction purposes were flat-fielded, and combined after re- 
moving the outliers for each pixel, which removes any stars 
present in the individual images. The result was subsequently 



' http://www.not.iac.es/instruments/notcam/calibration.html 
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Fig. 2. Light curves for our WFC r- and I-band observations. The best fit transit models are overplotted. The r-band data were 
detrended on a point-by-point basis (see text), and both hght curves have been corrected for airmass. The errorbars shown are set to 
the scatter in the out-of-transit basehne. 



subtracted from the science images, removing most of the struc- 
ture in the background. 

In the next step, aperture photometry was performed on 
GJ1214 and the reference stars, using an aperture of 10 pixels 
for the NOTCam dataset. The residual background was deter- 
mined by averaging the background signal in an annulus cen- 
tered on the stars with an inner radius of 30, and an outer radius 
of 50 pixels. Individual pixels in the annulus were clipped at 5-cr 
from the mean to prevent outliers from affecting the background 
measurements. 

The WHT-LIRIS Kc-band data were reduced in a similar 
way, but we constructed our own non-linearity curve from a se- 
quence of dome-flats at different exposure times. After the sky- 
subtraction, aperture photometry was performed with a 25 pixel 
radius, using an annulus from 30 to 50 pixels to determine the 
residual background, clipping the outliers in the same way as for 
the NOTCam data. 

The NOTCam Kj- and LIRIS K^-band light curves nor- 
malised by those from the reference stars are shown in the top- 
panels of Fig. |4] Further corrections for systematic effects are 
done simultaneously with the transit fitting and are discussed in 
Sect.lO 

4. Transit fitting 

4.1. Optical transits 

Since our goal is to measure the effective radius of the planet as 
a function of wavelength and not to determine the other transit 
parameters, we fit the transit at each band independently, keeping 
the imp act parameter and semi-major axis fixed to the values 
used by iBean et aTl (l2010l) . This has the added benefit that we 
can compare our results directly to previous literature values. To 



model the transit curve we use the formalism of iMandel & Agol 
(2002) with non-linear limb-darkening coeflicients as calculated 
by Claret (2000) for the INT I-band observations and by Claret 
(200|^r the observations at the other wavelengths. We used the 
limb-darkening parameters for the gridpoint closest to the stellar 
parameters from Charbonneau et al. (2009). The limb-darkening 
coeflicients are given in Table |2] We performed a Markov-Chain 
Monte Carlo (MCMC) analysis to find the planet-to-star radius 
ratio and the time of mid-transit, creating 5 chains with a length 
of 200,000 points for each band. From each of the chains the first 
20,000 points were trimmed to ensure that the initial conditions 
did not influence the r esults. After checking tha t the chains were 
well mixed using the i Gelman & RubinI ( Il992h test, we merged 
them. The best fit models are overplotted on the hght curves in 
Figs.|2]and[3] 

To assess the impact of red-noise on uncertainties of the fit- 
ted parameters, we used the residual permutation method, which 
uses the data itself to assess the uncertainties on the measured 
parameters due to correlated noise. The best fit model is added 
to the residuals from the MCMC fit, after the residuals were 
shifted by N points. The points are wrapped around when per- 
forming the shift. The new light curve is then refitted. F or the 
fitting we used the IDL MPFIT packag^ ( Markwardtll200 9). To 
assess the uncertainties on the parameters we use the range be- 
tween 16% and 84% of the distribution of the parameters for the 
+ l-cr uncertainty interval. In the INT r- and I-band observations, 
the uncertainty estimates from the residual permutation method 
are lower than those obtained from the MCMC analysis, this is 
probably due to the low number of frames that were obtained. 
In this case we have adopted the uncertainties from the MCMC 



" http://purl.com/net/mpfit 
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Fig. 3. The light curves for our GROND observations with the best fit models overplotted. All light curves have been corrected for 
a linear trend in the baseline. The errorbars shown are set to the scatter in the out-of-transit baseline. 



analysis rather than those from the residual permutation analy- 
sis. 

The fitted radius ratios for all bands (see also Table |3]l, are 
0.1143+0.0018 and 

0.1162+0.0005, for our INT r and I-band measurements, 
and 0.1198+^°j;26^ 0.1168+0.0010, 0.1162+0.0013 and 
0.1165+0.0013 for our GROND g-, r-, i- and z-band ob- 
servations respectively. The times of mid-transit are given in 
Table n 

4.2. Near-infrared transits 

Firstly, outliers were removed from the Ks-band light curve 
by excluding all points more than 1.2% away from a median 



smoothed light curve. In this way 25 points were removed. A 
clear residual trend in the out-of-transit baseline is visible in 
the light curve (Fig. |4]i, but we find no significant correlations 
with the position on the detector or airmass , which we often see 
in other near-infrared measurements (e.g. Ide Mooii & SnellenI 
2009; de Mooij et al. 2011). We therefore fitted the hght curve 
with a second order polynomial simultaneously with the tran- 
sit parameters. As for the optical light curves, we only fitted for 
the time of mid-transit and the planet-to-star size ratio, keeping 
the impact parameter and semi-major axis fixed to the values 
used by B ean et al, (2010) . The fits were again performed using 
an MCMC analysis, as for the optical data, using 5 chains of 
200,000 steps, trimming away the first 20,000 points and check- 
ing that they were well mixed CGelman & Rubin 1992) . Non- 
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Fig. 4. Light curves for our NOTCam Kj-band and LIRIS Kc-band observations. Top panels: Raw, normalised light curves with best 
fit model for the transit and systematic effects overplotted. Middle panels: Light curves coiTected for systematic effects, with the 
best fit transit model overplotted. Bottom panel: residuals after subtracting the best fit models. 



linear limb-darkening parameters from lClarefl (l2000h were used 
(given in Table |2]i. The best fit model is overplotted in the top- 
left panel of Fig. |4] The light curve with the corrected baseline 
is shown in the left middle panel, and the residuals in the bottom 
left panel of the same figure. 

Just as for the optical data, the impact of correlated noise 
is investigated using the residual permutation method. The es- 
timates for the uncertainties given by the residual permuta- 
tion method are almost identical to those we obtain from our 
MCMC analysis, which indicates that, apart from the trend in 
the baseline, the systematic effects are small compared to the 
random noise. We measure a radius ratio in the Kj-band of 
0.1189+0.0015. 

The LIRIS K^-band light curves were fitted in a similar way 
to the NOTCam Kj-band data. First, 13 outlier points were re- 
moved from the light curve. In addition to the systematic trends 
seen in the baseline, which again do not correlate with airmass 
or with position on the detector, we found a coiTelation with 
the fraction of the total flux of the star contained in the bright- 
est 1 1 pixels, which we take to be a proxy for the seeing (note 
that due to the strong defocus we cannot measure the seeing di- 
rectly). We therefore fitted the relevant transit parameters simul- 



taneously with a second order polynomial and this seeing proxy. 
We used a MCMC analysis similar to that used for the Kj-band 
light curve. In Fig. |4] the best fit model is overplotted in the top 
right panel, and the light curve coiTected for systematic effects 
is shown in the middle right panel. The residuals from the best 
fit model are shown in the bottom right panel of the same figure. 

The estimates for the uncertainties given by the residual per- 
mutation method are ~50% higher than those that were obtained 
from the MCMC analysis, and there is also an asymmetry be- 
tween the upper and lower uncertainties. 

As can be seen in the bottom right panel of Fig. |4l a sig- 
nificant feature is visible in the coiTected K^-band light curve 
between phase -0.012 and phase -0.006. Although this feature 
could be due to starspots, at the same moment a sudden increase 
in the difference in background occurs between the target and 
the reference star. At the moment of this 'bump' we also see a 
change in the shape of the PSF at lower intensities across the 
field of view of the detector (especially around the diffraction 
spikes caused by the support of the secondary miiTor, see Fig.|5]l. 
It seems that instrument-related effects influence the measure- 
ments at this point. To check the influence of this feature on 
the fit to the light curve, we excluded it and refitted the time- 
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Fig. 5. Images of the PSFs for our LIRIS Kc-band observations 
taken at two different times. The left column is for GJ1214, while 
the middle and right columns are for the two reference stars. The 
top panels show the average PSF in a sequence of 3 1 frames just 
before an anomalous feature is visible in the light curve, while 
the middle row shows the average PSF during the anomalous 
feature. The difference between the PSFs is shown in the bottom 
row. 



series. While for the light curve as a whole we find a radius 
ratio of 0.1162+° °°'^, we find a radius-ratio of 0.1187±0.0018 
(~1.5cr deeper) when this feature is excluded from the fit. Since 
we do not want to bias our results, we use the radius-ratio ob- 
tained from fitting the entire K^-band light curve for the rest of 
the paper The extra uncertainty is accounted for by increasing 
the error on the radius ratio by 50%. 

5. Stellar variability 

ICharbonneau et all (l2009 l) found that the star GJ1214 is variable 
at the level of 1% in the MEarth bandpass, with an apparent pe- 
riod of ~80 days, attributed to the rotation-modulated variations 
in the star-spot levels. Recently Berta et al. (2011) presented a 
stellar variability analysis covering several years of data, find- 
ing a period of 52.3+5.3 days. This star-spot induced variability 
can have a major impact on transit photometr y taken at different 
epochs and/or at different wavelengths (e.g. ICarter et al.ir201 U 
ISing et ani201 Ibh . Since our data are taken over the course of 
3.5 months, and our aim is to measure differences in the transit- 
depth as a function of wavelength, both the temporal and spec- 
tral effects from starspots are present and need to be taken into 
account. 

If a star-spot is present on the stellar surface during transit, 
firstly the planet could transit the spot. Because the star-spots 
have a lower surface brightness than the unspotted stellar sur- 
face, a spot that is occulted by the planet would make the transit 
less deep and bias the planet-to-star radius ratio to lower values. 
In the case the planet does not cross the star spot, the transit ap- 
pears deeper, since the average surface brightness along the path 
of the planet is higher than that of the spotted surface, leading 



to a decrease of the effective radius of the star. Unocculted spots 
will therefore bias the planet-to-star size ratio to higher values. 
In addition, since the star spots have a lower effective tempera- 
ture than the stellar surface, and are therefore relatively red, these 
effects will be stronger in the blue part of the optical spectrum 
than in the near-infrared. An expr ession for the bias due to un- 
occulted starspots can be found in lSing et al.l (1201 Ibl) . for more 
generalised equation s for these corrections we refer the reader to 
Desert et all (1201 Ibh . 

5. 1 . Correcting for the stellar variability 

Changes in the spot fraction can be monitored by measuring the 
stellar brightness as a function of time, which can subsequently 
be used to link the observations taken at different epochs. Note 
that this does not take into account the base-level of spots, which 
is unknown but can still significantly influence the measured 
planet-to-star radius ratio as a function of wavelength. This is 
discussed in Sect. 16.41 

To link our INT I-band observations to our INT r-band ob- 
servations, we obtained a set of r-band out-of-transit frames on 
July 29, 2010, in order to measure the change in stellar bright- 
ness between the two nights. These r-band frames were reduced 
in the same way as the transit observations, and aperture pho- 
tometry was performed on the same set of stars as for the ob- 
servations on May 26. We find that between May 26 and July 
29 the star has decreased its brightness by AF=2%. We therefore 
correct our INT r-band and NOTCam K^-band points to match 
the stellar brightness of the INT I-band observations. We use the 
expression for the effects of starspots on the planet-to-star ra- 
dius rati o from Sing et al. ( 201 lb) , where we use the NextGen 
models (iHauschildt et alJll999h for flux of both the star and the 
spots. We assume a spot temperature of Tjpo,=2800 K, which is 
200 K cooler than the effective temperature of the star We find a 
correction ARp/R, =0.001 1 in the r-band and AR^/R, =0.0003 in 
the K 5 -band. 

We linked the stellar brightness during the night of the 
GROND observations to that during the INT r-band observa- 
tions using the data from Berta et al. (2011). Subsequently, this 
was linked to our INT I-band observations using the offset found 
above. These two steps cancel each other out within the un- 
certainties and therefore no correction was performed on the 
GROND measurements. We also did not correct the K^-band 
observation, since we have no flux measurements of the stellar 
brightness in the r-band on that night. However, the correction is 
expected to be small (AR^/R* <0.0003), and will therefore not 
significantly influence the measured transmission spectrum. The 
corrections for the radius ratios are given in Table [3j 

6. Discussion 

6.1. The transmission spectrum of GJ1214b 

Using the measured planet-to-star radius ratios in all of our 
bands, corrected for stellar variability, we can construct a trans- 
mission spectrum for GJ1214b, which is shown in Fig. |6l All 
data points are consistent with a constant planet-to-star size ratio 
within Icr, except for the g-band point at 0.46/im, which deviates 
by ~2cr. 



6.2. Atmospheric models 

To investigate what constraints our observations can place on the 
atmosphere of GJ1214b, we compare our results with theoretical 
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Table 1. Parameters used for the transit fitting of GJ 12 14 and the measured times of mid- transit. 



Parameter 




Value 


Reference 


a/R. 




14.9749 


(1) 


b=a/R,cos(i) 




0.27729 


(1) 


P 


(days) 


1.580408346 


(1) 


T 

A star 


(K) 


3030 


(2) 


Tc -2,400,000 (HJD) 


(INT r) 


55342.66081 ±0.00029 


(3) 




(NOT K,,) 


55342.66057±0.00023 


(3) 




(GROND g) 


55380.59049±0.00040 


(3) 




(GROND r) 


55380.59045±0.00019 


(3) 




(GROND i) 


55380.59044±0.00022 


(3) 




(GROND z) 


55380.59013±0.00017 


(3) 




(INT I) 


55407.45754±0.00007 


(3) 




(WHT K^) 


55426.42274±0.00032 


(3) 



References. (D lBean et aLlfcOlOh : (2) ICharbonneau et alJ l l2009l) : (3) this work. 

Table 2. Limb-darkening coefficients used for the fitting of the transit light curves. 



Instrument 


Band 


Ci 


Ci 


Ci 


Ci 


Reference 


GROND 


g 


0.2306 


1.0639 


-0.3180 


0.0119 


(2) 


GROND 


r 


0.5654 


0.1274 


0.5991 


-0.3119 


(2) 


WFC 


r 


0.5654 


0.1274 


0.5991 


-0.3119 


(2) 


GROND 


i 


0.6819 


0.4526 


-0.1692 


-0.0041 


(2) 


WFC 


I 


0.6061 


0.8676 


-0.7015 


0.1881 


(1) 


GROND 


z 


0.9613 


0.5228 


-0.8702 


0.3438 


(2) 


NOTCam 


K, 


1.9371 


-2.8039 


2.3942 


-0.8098 


(1) 


LIRIS 


K, 


1.9371 


-2.8039 


2.3942 


-0.8098 


(1) 



References. (n lclarea ( l2000l) : (2) Ictoea ( 120041) . 
Table 3. Radius ratios from our observations 



Instrument 


Band 


A, 
(jum) 


V / measured 






GROND 


g 


0.46 


1198+0"026 




-0.0007 


GROND 


r 


0.62 


0.1168±0.0010 




-0.0005 


WFC 


r 


0.63 


0.1143±0.0018 


0.0011 


-0.0006 


GROND 


i 


0.76 


0.1162±0.0013 




0.0000 


WFC 


I 


0.82 


0.1162±0.0005 




+0.0005 


GROND 


z 


0.90 


0.1165±0.0013 




+0.0009 


NOTCam 


K, 


2.15 


0.1189+0.0015 


0.0003 


+0.0023 


LIRIS 


K, 


2.27 


0.1162+0.0030 




+0.0021 



Notes. Radius ratios as determined from our observations in dilferent filters (column 4), the correction that is applied to the radius ratios due 
to stellar variability (column 5), as well as the wavelength-depedent change in radius ratio assuming a background of unocculted starspots that 
cover 10% of the stellar surface (column 6) (note that the difference in i-band is set to 0) (see also Sect. 16.41 ). The instruments, filters and central 
wavelengths are indicated in columns 1 through 3. 



models of the planet's atmosphere. The models include Rayleigh 
scattering as well as molecular absorption features and collision- 
induced absorption from H2-H2 (Borvsow 200 2). For the trans- 
mission spectrum we calculated the light passing through layers 
of atmosphere between 110 '' and 5 bar At higher pressures we 
consider the atmosphere to be fully opaque. At each point we as- 



sume that the atmosphere can be described by an average profile 
in hydrostatic equilibrium. For the temperature-pressure profile 
we use one that is similar to that used bv Miller-Ricci & FortnevI 
( l2010h . Our models include absorption b y water, for which the 
opacity data from the HITEMP database dRothman et alJl2010l) 
were used. In addition we also include methane, for which 



8 



De Mooij et al.: Transit observations of GJ1214b 



0.135 



0.130 - 



0.125 - 



0.120 



0.115- 



0.1 10 



0.5 



i I 
i > 




\ I 

■\ ' / _ 



1.0 1.5 

Lombdo (^m) 



2.0 



2.5 



Fig. 6. Transmission spectrum of GJ1214b from our data. Overplotted are models for the atmosphere of GJ1214b. The green 
(dashed-dotted) line is a hydrogen rich atmosphere with a solar metallicity, the red (dashed) curve is for an atmosphere that is 
hydrogen rich with a sub-solar metallicity and a cloud layer at 0.5 bar, the blue (dashed-triple dotted) line is for an atmosphere that 
is dominated by water. The dotted curves at the bottom of the plot show the various transmission curves of the filters used for our 
measurements. 



the data from HITRAN 2008 ("Roth man et"al]|2009l) were used. 
Water and methane are expected to be the main absorbing gases 
in the atmosphere. A line-by-line code was used to calculate the 
opacities for the different gases, assuming a homogeneous mix- 
ing of all the different species present in the atmosphere. We 
used a Voigt profile for the individual lines and applied a line 
wing cut -of f of 30 cm~' to simulate a sub- Lorentzian line pro- 
file (e.g. lBaiiev & Kedziora-Chudczed201 ih . This gives rise to a 
lower continuum and more pronounced absorption features than 
GJ1214b model results presented elsewhere. 

We have chosen to generate three models to make a qualita- 
tive study of the atmosphere of GJ1214b, and we have overplot- 
ted these models on the observed transmission spectrum shown 
in Fig. |6l The models were matched to the mean of the mea- 
sured radius ratios between 0.7 fim and 1 fj.m. The first model, 
the green (dash-dotted) line in Fig. |6l is for a (geometrically) 
thick atmosphere with a solar composition and no cloud layers. 
The concentrations of water and methane are 310 The dom- 
inant hydrogen gives rise to a high atmospheric scale height, as 
can be seen by the strong molecular features in the near infrared. 
Fig.|6]indicates that a solar metallicity atmosphere gives features 
that are too strong compared to our measurements, resulting in a 
of 26.7. The second model, the red (dashed) line in the figure, 
is for an atmosphere with a sub-solar metallicity, and includes 
a grey cloud layer at a pressure of 0.5 bar. The concentrations 
of water and methane are 5.6-10"^ and 5-10"^ respectively. This 
model is consistent with our measurements (j(;^-3.2). The third 
model, the blue (dashed-triple dotted) line in the figure is for a 
water dominated atmosphere. This model has its molecular fea- 
tures significantly suppressed due to the decreased scale height. 
This model is also consistent with all our observations, but with 
a higher of 15.1 compared to that of the second model. This 
is mainly due to the g-band measurement. The latter two mod- 
els are chosen such that they can reproduce the relatively flat 
transmission profile seen in the red part of the optical spectrum. 



These values can be compared to that of a straight line, which 
has ;t'2^ 10.5. 

The data can be well fitted with a model with a sub-solar 
metallicity and a cloud layer. The major disagreement between 
our data and the model for a water-dominated atmosphere comes 
from the GROND g-band measurements, which could be due to 
Rayleigh scattering, and to a lesser extent from the NOTCam 
Kj-band observations. It is therefore important to obtain further 
observations to confirm these measurements. 

6.3. Comparison with previous measurements 

During the past year there have been m any transit measure ments 
of GJ1214b at different wavelengths. iBean et aP (1201 Ol) mea- 
sured the transit-depth in 11 wavelength bins from 790 nm to 
1000 nm. Our radius ratios obtained in the same wavelength 
range are consistent with their measurements. iBeanetal J (120101) 
cannot explain their observations with a thick, hydrogen-rich at- 
mosphere, and therefore favour models for a water-dominated 
atmosphere of GJ1214b. However, they cannot rule out clouds as 
an explanation for their relatively featureless transmission spec- 
trum. This is also seen when comparing their data with our mod- 
els, the x^ for the water model is 9.1, while that for the model 
with a sub-solar composition with a cloud layer is 7.8, and for a 
str aight line x^ = 6.0. 

ICroU et aP (|201 ih have presented measurements of the 
transit-depths in the J- and K^bands, as well as in a band cen- 
tered on a methane feature at 1 .69 fim. Their measurements in 
the J- and K^-band provide evidence for an increase in the plan- 
etary radius at 2.15 //m compared to that at 1.15 fim. This in- 
crease is consistent with the presence of a thick, hydrogen-rich 
atmosphere. The additional measurement of the transit-depth 
at 1.69 fim, which probes a methane feature in the H-band, is 
lower than that observed in the J-band. This point however has 
a large uncertainty, and is therefore still consistent with an at- 
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mosphere with a low mean-molecular weight, but only if the 
methane abundance is low. Comparison of these data with our 
models shows a significant increase in the between the sub- 
solar metalli city model and th e water model of A;t'^=20.1, as 
advocated bv lCroU et alj (1201 lb . The large depth in the K^-band 
for the sub-solar metallicity model is due to H2-H2 CIA and not 
due to methane. 

In con trast, using t h e IRAC camera on the Spitzer Space 
Telescope, iDesert et alj (1201 lab find radius ratios in the IRAC 
ba nds at 3.6 and 4.5 u m that are consistent with the observations 
of iBean et alJ (1201 Oh . showing a flat transmission s pectrum ex- 
tend ing into the mid-infrared. The observations by Dese rt et"al] 
(l2011a) therefore reject a cloud-free methane rich atmosphere 
because it would require a significantly larger radius at 3.6//m 
th an is observ e d. With in the Icr errorbars, the observations 
of iDesert et al.l (1201 la l) are consistent with both the sub-solar 
metallicity model with a large scale-height and with the model 
for a water-dominated atmosphere. 

In addition to ou r own g-band measurement, observations 
bv 'Carte r et all (1201 ih in the r- and z-band show evidence for an 
increase in the radius-ratio of the planet when going from the 
z-band to the r-band at the Icr level. However, they also caution 
that stellar spots are important when comparing observations at 
different dates and wavelengths. The measured increase of the 
planet-to-star radius ratio is 0.00184+0.00166, which is con- 
sistent with Rayleigh scattering, where the expected difference 
is ~0.0015. It is, however, difficult to dire ctly compare these 
measurements with our measurements, since^ Carter et alj(l201 Ih 
used a different impact parameter and semi-major axis. 

In the top panel of Figs. [T] and [8] we show all the cur- 
rently available measurements of the transmission spectrum of 
GJ1214b that were obtained using the same impact parameter 
and scaled semi-major axis. We caution, however, that when 
making a direct comparison between literature values and our 
observations one has to be careful, since these were taken at dif- 
ferent epochs, and therefore at different levels of stellar variabil- 
ity. 

Ignoring these subtle effects, the available data favour the 
model for a hydrogen-rich, low metallicity atmosphere with a 
cloud layer at 0.5 bar {x^-\A3) over a water dominated at- 
mosphere (x^-A^.l) and a flat line (x^-A6.%). This is mainly 
due to our g-band measurements and t he measurements in the 
Kj-band presented in this paper and in 'Croll et all (1201 i h. The 
Spitzer measurements from .Desert et al., (,201 lal) are consistent 
with both models. 

6.4. The impact of unocculted starspots 

As discussed in Section |5] unocculted spots can also produce a 
wavelength dependent signal in the transmission spectrum, by 
altering the eff'ective radius of the star. Since only the variation 
in the starspots is traced by the variability measurements, and 
not the total level of spots, we have no constraints on the base 
level of spots present on the stellar surface. There is a priori no 
reason to assume that the only spots present on the stellar surface 
are those responsible for the variability. This would result in flat 
regions in the light curve when no spots are visible on the stellar 
hemisph ere pointed t oward Earth, which is not seen in the light 
curve bv lBerta et all (1201 lb . 

From this variability curve we know that GJ1214 varies in 
brightness at the ~l% level in the MEarth bandpass. This corre- 
sponds to a variation in the spot-fraction of ~2.5%, assuming a 
spot temperature of 2800K. We might expect that this variability 
is the result of a distribution of many small spots on the stellar 



surface. This is sup ported by th e fact that for 2 out of the 16 
transits observed bv lCarter et alJ (1201 1 *) the planet crosses over 
a small star spot. Since these spots are possibly distributed in a 
random pattern in longitude on the stellar surface, we argue that 
the fraction of the star that is covered in spots is significantly 
higher than the 2.5% necessary to explain the variability. 

To investigate at what level the observed transmission spec- 
trum can be influenced by a base level of unocculted spots, 
we calculated the wavelength dependent bias in the radius ra- 
tio for several assumed base levels, which is overplotted on 
the observed transmission spectrum in the middle panels of 
Figs. Q and |8| This bias was calculated in the same way as the 
coiTections for stellar variability in Sect. 15.11 However, instead 
of the observed stellar variability, we calculated the bias for dif- 
ferent amounts of stellar dimming in the GROND i-band. As can 
be seen it is very unlikely that star spots are responsible for the 
higher g-band and Kj-band values. To explain the g-band mea- 
surement, a large part of the stellar surface (>20%) has to be 
covered in starspots, which would lead to a much lower radius 
ratio in the near and mid-infrared than observed. 

However, this argument can also be turned around. Although 
the observed transmission spectrum cannot be explained by a 
baseline of unocculted starspots, if there is a significant baseline 
present, the shape of the planet's transmission spectrum can also 
be altered in such a way that the observed transmission spec- 
trum appears relatively flat. We therefore also investigate how 
the starspots can change the shape of the transmission spectrum. 
To demonstrate this impact, we applied a correction for the spot 
base level, assuming that the star is dimmed by 5% in the i-band, 
which corresponds to a spot covering fraction of 10% of the 
stellar surface. In the bottom panels of Figs. |2|and[8]we show 
this corrected transmission spectrum, where, to make compari- 
son with the earlier results easier, we have shifted the values such 
that the i-band radius ratio is unaltered. The effect of this correc- 
tion is to push the near and mid-infrared points to larger radius 
ratios with respect to the optical measurements. The increased 
radius ratio in the infrared is then more consistent with a geomet- 
rically thick atmosphere, although we require a strongly reduced 
methane content in order to explain both our K^-band and IRAC 
points. At the red part of the optical spectrum the correction 
introduces a tilt towards smaller radii for shorter wavelenghts, 
which is opposite to the tilt expected from Rayleigh scattering 
in the atmosphere of the planet, which causes an increase in the 
radius ratio. 

From this example it is clear that not only stellar variability 
needs to be taken into account when comparing the transmission 
spectrum at different wavelenghts and times, but that it is also 
extremely important to investigate the impact of the base level 
of star spots on the observed transmission spectrum. 

7. Conclusions 

We have presented observations of 8 transits of GJ1214b, in the 
g-, r-, i-, I- and z-bands in the optical, and K,- and K^-bands in 
the near-infrared, allowing us to determine the planet-to-star ra- 
dius ratio of this system at these wavelengths. All measurements 
show within errors a similar size ratio, except for the g-band, for 
which this ratio is larger at the ~2cr level. If real, this could be 
attributed to Rayleigh scattering within the planet's atmosphere, 
and subsequently points to a large scale height and therefore a 
hydrogen/helium dominated atmosphere. However, this requires 
more observations to confirm. 

When combining our observations with all the currently 
available literature values, the data again favour a model with an 
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extended atmosphere, a low methane content and a cloud-layer 
at 0.5 bar, above that for a water dominated atmosphere. 

In addition, it is very important to take into account the base- 
level of stellar spots, which can significantly alter the shape of 
the transmission spectrum, by increasing the radius ratios mea- 
sured at optical wavelengths compared to those at near- and mid- 
infrared wavelengths. When this effect is removed, the relative 
increases in the planet-to-star size ratios in the near- and mid- 
infrared make those measurements more consistent with an at- 
mosphere with a large scale-height compared to an atmosphere 
which is dominated by water vapour. 
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Fig. 7. Top panel: Transmission spectrum of GJ1214b. Overplotted are the atmospheric models for GJ1214, using the same linestyles 
as in Fig. |6]for the atmospheri c model s. Our data are plotted as fill ed circles, the VLT data from Bean et al. (2010) as crosses, the 
Spitzer data from Desert et al.l (l201 lal) as stars, and the data from ICroll et al.l (|201 ih as diamonds. Middle panel: The same data 
points as in the top-panel but now with different models overplotted showing the change in radius ratio due to different amounts 
for the base levels of stellar spots in the GROND i-band. The dotted, dashed, dash-dotted and the dashed-triple dotted lines are for 
10%, 5%, 2.5% and 1% contribution from stellar spots in the i-band. Bottom panel: The observed transmission spectrum, but now 
corrected for a spot-dimming in i-band due to spots of 5% (The dashed curve in the middle panel), which corresponds to a spot 
coverage of 10% of the stellar surface. The i-band radius ratio was kept at the same level as in the other panels. 
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Fig, 8. Same as Fig. |7]but now only showing the optical part of the transmission spectrum 



0.8 



0.9 



1.0 



13 



